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Enhancement of Antiferromagnetic Correlations below Superconducting Transition 

Temperature in Bilayer Superconductors 
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Motivated by the recent experiment in niultilayered cuprate superconductors reporting the en- 
hancement of antiferromagnetic order below the superconducting transition temperature, we study 
the proximity effect of the antiferromagnetic correlation in a bilayer system and also examine the 
possibility of a coexistence of antiferromagnetic order and superconductivity. We present the result 
of mean field theory that is consistent with the experiment and supports the proximity effect picture. 

PACS numbers: 74.78.Fk, 74.72.-h 



I. INTRODUCTION 



In the study of high-temperature superconductivity in 
the cuprates, there has been intense interest in the inter- 
play between antiferromagnetism and superconductivity. 
The undoped parent compound, which is a Mott insu- 
lator, is an antiferromagnetic (AF) long-range ordered 
state. Upon carrier doping, the AF state is converted 
into the high-temperature superconducting (SC) state. 
In the mechanism of superconductivity, the AF correla- 
tion is believed to play an important role. In order to 
investigate the AF correlation effect on superconductiv- 
ity, multilayer systems are useful. In multilayer systems, 
the number of Cu02 layers per unit cell, n, takes n > 2. 
From systematic studies of the n dependence of SC tran- 
sition temperature Tc, it is found that the maximum Tc 
is obtainedi— for the case of n = 3,4. Since coupling 
between the Cu02 planes within the unit cell is stronger 
than that between the Cu02 planes in different unit cells, 
the AF correlation can be strong in multilayer systems. 
There is also charge imbalance over the Cu02 layers as re- 
vealed by the ^^Cu Knight shifts, and thus the multilayer 
systems form a natural heterostructure at an intermedi- 
ate carrier-concentration regioui^ 

From nuclear magnetic resonance (NMR) experiments 
it was suggested that AF moments survived even in a 
SC phase4ri A possibility of a phase separation is ruled 
out because the NMR signal associated with the para- 
magnetic state was not observed)^ An important issue 
here is whether SC order coexists with AF orde r^i^°i^^ 
or not. Recently, Shimizu et al. studied the temper- 
ature dependence of the AF moment in a five-layered 
cuprate Ba2Ca4Cu50io(F,0)2 by NMR measurements J^ 
They reported that the AF moment in the outer layer is 
enhanced below Tc. Such an enhancement of the AF mo- 
ment below Tc docs not appear in the coexistence phase 
of other superconductors such as iron arsenide supercon- 
ductor Ba2(Fei_2;Co2;)2As2, where the AF moment re- 
duces at Tcji^ Therefore, the enhancement seen in the 
five-layered cuprate is not necessarily common to all of 
the coexistence phase of superconductors. Noting that 
the iron arsenide is a single-layer system, we may spec- 
ulate that the multilayer nature is essential for the en- 



hancement, and thus need to clarify the mechanism of the 
AF moment enhancement triggered by the SC transition 
under the presence of multilayers. 

In this paper, we study the interplay between SC and 
AF orders within a mean field theory. To describe a 
multilayer system, where there arc two types of layers 
associated with outer planes and inner planes, we con- 
sider a bilayer system with different electronic correla- 
tions in each layer. By controlling interaction parameters 
for the AF correlation and the SC correlation, we study 
the proximity effect between two layers with different or- 
der and also examine the possibility of coexistence of SC 
and AF orders. We find that the proximity effect leads 
to an enhancement of AF order below the SC transition 
temperature that is consistent with the experiments^ By 
contrast, we find qualitatively different behaviors of the 
order parameters for the coexistence case. 

This paper is organized as follows. In Sec|lll we in- 
troduce the bilayer Hamiltonian and describe the mean 
field theory. In our model, a coexistence phase of SC 
and AF orders is possible within a single layer. This 
coexistence phase is described in Sec lIIIl In Sec lIVl we 
describe the results about the proximity effect in the bi- 
layer system. We show that the temperature dependence 
of the AF moment is consistent with the experiments^ 
We also examine coexistence phases and show that the 
temperature dependence of the order parameters is quite 
different from the experiment. Finally, Sec|V] is devoted 
to summary and discussions. 



II. MODEL AND FORMALISM 

We consider a bilayer system with interactions which 
stabilize SC and/or AF order. The two layers are cou- 
pled through an interlayer tunneling. The Hamiltonian 
is given by 



n=J2ni + n±, 



1=1.2 



where the interlayer tunneling term Ti^ is 

4,fe^C2,fc^+H.C, 



'H± = -ip X! ' ^' 



(1) 



(2) 



k,a 



Here c[^^ {ci^ka) creates (annihilates) electrons with in- 
plane momentum k and spin a at layer /. We assume 
that the interlayer hopping matrix tp is independent of 
k. The Hamiltonian for the I layer is 



fc,CT j 

-91 Yl /('^)/(^')4fcT4-fciC(,-fcuc;,fc't, (3) 

with ^fc = — 2i(cosfca; + cosky) — fj,. (Hereafter we set 
the lattice constant to unity.) The hopping of electrons 
within each layer is restricted to the nearest neighbors 
given by t. Energies are measured in units of t in the 
following analysis. We focus on the half-filing case in 
each layer so that we set the chemical potential fi ~ 0. 
This means that we neglect charge redistribution between 
the two layers. The second term of Hi describes the AF 
interaction. The operator cj- (cija) creates (annihilates) 
electrons at site j. The third term of Hi with f{k) = 
{coskx — cosky) / 2 describes the interaction for dx^^^y^- 
wave pairing. 

We define the d-wave SC order parameter in the / layer 
with N sites by 



^'-^E^'W' 



(4) 



with 



Ai(fc) = f{k){ci,k^cu-ki). (5) 

The AF order parameter in the I layer is defined by 



mi 



27V 



fceRBZ 



where the summation with respect to k is taken over the 
reduced Brillouin zone lA:^:! -I- |A:j,| < tt and the nesting 
vector is <5 = (tt, tt). By using the order parameters, the 
mean field Hamiltonian at the I layer reads 

'^f = Y. Cl^MikCik + ^NVim^ + Ngi\^i\', (7) 

fceRBZ 



where Cik = ( Q,fct Q,fc+Qt c^^ki 4,-k-Qi ' and 



A//fc = 



?fc -^miVi -giMk) 

-4m, V/ ik+Q -geAi{k) 

uMI^T -a -"^miVi 

-geAiik)* -AmiVi -^k+Q 



(8) 
Based on the mean field Hamiltonian of the whole system, 
J2i ^r^ + H±,we solve the mean-field Eqs. (g]) and ^ 
numerically using the 100 x 100 discretized Brillouin zone. 



III. COEXISTENCE PHASE IN SINGLE LAYER 
SYSTEM 

Before going into the analysis of the bilayer system, 
we examine a coexisting phase in a single layer system 
described by Hf^^. Reflecting the difference in symmetry 
between the AF gap created by to/ 7^ and the SC gap 
A/, the coexistence phase of SC and AF orders can be 
stabilized»i^ii^ The situation is similar to the slave-boson 
mean field theory of the t-J modeU^ and the Hubbard 
model in the strong-coupling limit.— 

Figure \V[^a,) shows the parameter range oi g = gi for 
the coexistence phase at T^ = 14 = 0.5. For g < 5.0 the 
system is a pure AF state while for g > 6.3 the system 
is a pure SC state.— So the coexistence phase appears for 
5.0 < g < 6.3. For the case oi V = 0.4, the parameter 
for the coexistence changes as 3.7 < g < 4.5. In order to 
confirm that the state with A 7^ and m 7^ is the global 
minimum of the free energy, we computed the following 
energy at T = for < A < 0.2 and < m < 0.25: 



E = 



E ^" 



+ 4:NVm^ +Ng\Af' 



(9) 



Q,B„<0 



where a runs over the all eigenstatcs of the mean field 
Hamiltonian and Ea are the eigenenergies. We examined 
several cases and confirmed that the coexistence phase 
solution corresponds to the global minimum of the en- 
ergy. We examined the s-wave case as well, but there is 
no coexistence phase. 

There are two types of coexistence phases. One is the 
phase with strong AF order and weak SC order, resulting 
in the AF transition temperature Taf > ^c, and the 
other is the phase with Taf < 7c. Figure [Ijb) shows the 
temperature dependence of the order parameters in the 
coexistence phase with Taf > Tc at g = 5 and V = 0.5. 
For comparison, the pure AF case at 5 = and V — 0.5 
and the pure SC case at g = 5 and V^ = are also shown. 
For this choice of the parameters. Tap is slightly higher 
than Tc- Therefore, the system first exhibits AF order 
upon decreasing temperature. T^ is somewhat reduced 
because of the presence of AF order. The occurrence of 
the SC order also affects the AF order: the temperature 
dependence of m in Fig. [Ijb) deviates from the pure AF 
case at Tc, resulting in the reduction of m. This behavior 
is in contrast to the case of the multilayer cuprates where 
the enhancement of the AF order is observed J^ We note 
that, in the coexistence phase with Tc > Tap, Tc is the 
same as the value of the pure SC case but Taf is reduced 
from the pure AF case value. 



IV. PROXIMITY EFFECT IN BILAYER 
SYSTEM 

In this section, we study the bilayer system. Our pur- 
pose is twofold. First, we study the proximity effect in 
the bilayer system. Second, we examine the stability of 




FIG. 1. (Color online) (a) The AF order parameter m and 
the d-wave SC order parameter A versus the SC interaction 
parameter g at V = 0.5 in the single layer system at T = 0. 
(b) The temperature dependence of the order parameters in 
the coexistence phase for g = 5 and V = 0.5. For comparison 
we show the pure SC case (g = 5 and V — 0) and the pure 
AF case [g — and V = 0.5) as well. 



the coexistence phase in a single layer under the presence 
of intcr-laycr tunneling. 

To start with, wc examine the interlayer tunneling ef- 
fect. Depending on the value of tp, there are a strong 
tp regime and a wesik tp regime. Figure [2] shows the tp 
dependence of the AF order parameters mi and TO2 at 
Vi = 0.4 and V2 = 0.5 with g^ = g2 = 0. For tp < 0.22, 
we see that the values of the order parameters are not 
so much affected by the increase of tp. This weak tp 
regime is not suitable for describing multilayer systems 
because each layer is almost independent. In fact, the 
change of mi and 771,2 in the weak tp regime is described 
by the second-order perturbation theory with respect to 



tp. At tp = 0.22, there is a first order transition between 
the weak tp regime and the strong tp regime as shown in 
Fig. [21 For tp > 0.22, the order parameters exhibit strong 
tp dependence. In this strong tp regime, tp is larger than 
the excitation gap created by AF order. Therefore, the 
order parameters are reduced due to the change of the 
Fermi-surface topology. In the large tp limit, the nonin- 
teracting single-body electron states are well described by 
the bonding state and the anti-bonding state. The Fermi 
surface splits into two pockets centered at the F point 
and M point. Qualitatively similar behaviors are found 
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FIG. 2. (Color online) The AF order parameters ttii and 
7772 versus the interlayer hopping parameter tp in the bilayer 
system for V\ — 0.4 and V2 ~ 0.5. The inset is a schematic 
view of the system. 



in the tp dependence of the SC order parameters. In the 
following analysis, we focus on this strong tp regime and 
set tp ~ 0.3. 

Now we investigate the bilayer system. We consider 
three cases. In all cases, we assume 1^ = 0.5 and (72 = 
for / = 2. Therefore, the intrinsic order in the I = 2 layer 
is restricted to AF order. For the / = 1 layer, wc assume 
gi > 3 and Vi > 0. 

The three cases are (i) Vi = 0, (ii) Vi = 0.4(< V2), and 
(iii) Vi = 0.6(> ¥2). It turns out that the ground state 
of case (i) consists of SC order in the I = 1 layer and AF 
order in the I = 2 layer when gi > 3.7. The ground states 
of cases (ii) and (iii) are the coexistence state of SC and 
AF orders in the I ~ 1 layer and only AF order in the 
/ = 2 layer. Cases (ii) and (iii) are distinguished by the 
strength of the AF order in the two layers: 7^1 < 7^2 in 
case (ii), while m,i > m2 in case (iii). 

Figure [3lja) is a schematic view of case (i). The tem- 
perature dependence of the order parameters is shown in 
Figs.[3l^b)-(d) for different values oi gi. For gi = 3, there 
is no SC order in the I = 1 layer. Intrinsic AF order 7712 
appears in the / = 2 layer for T < Taf = 0.14 as shown 
in Fig. [3l^d). As a result of the proximity effect, a finite 
value of 7^1 is induced as shown in Fig.[3]Jc). This induced 
7771 decreases with increasing gi, because the intrinsic SC 
order by gi competes with the induced AF order. On the 
other hand, the intrinsic AF order 7^2 in the I = 2 layer 
increases with increasing gi. At gi ~ 4.0 there is a SC 
transition in the I = 1 layer as shown in Fig.[3Jb). In the 
presence of the non-zero SC order parameter, Ai, in the 
I ~ 1 layer, a finite value of A2 is induced in the / = 2 
layer (not shown) because of the proximity effect. The 
SC transition affects 7712. As shown in Fig. El^d), there 
is a clear enhancement of 77^2 below the SC transition 
temperature Tc = 0.06 for the case of gi = 4. A similar 
enhancement is also found in the case oi gi = 5 below 
Tc = 0.12. This enhancement of 7772 below T^ is consis- 
tent with the experiment in the multilayer cuprate J^ For 
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FIG. 3. (Color online) (a) Schematic view of the system for 
the case (i) (see the text). The temperature dependence of the 
SC order parameter Ai in layer-1 (b), the AF order parameter 
nil in layer-1 (c), and the AF order parameter 7712 in layer-2 
(d) for various gi. The vertical dashed line represents Tc in 
the case of gi — 4.0. 



(71 = 6 this enhancement is masked because Tc is higher 
than Tap- 

The enhancement of AF order below Tc is also found 
in the case (ii) (Fig. UJ^a)). Figures |3l^b)-(d) show the 
temperature dependence of the order parameters for dif- 
ferent values of gi. For the case of gi < 5.1, there is 
no intrinsic SC order in the I = 1 layer. For the case of 
gi = 5.2, there is a SC transition at Tc = 0.13. Below 
Tc, TO2 is clearly enhanced, although the enhancement is 
much reduced compared to the case (i). Similar behav- 
iors are observed for gi > 5.1. Again this behavior is 
consistent with the experimentji^ 



FIG. 4. (Color online) (a) Schematic view of the system 
for the case (ii) (see the text). The temperature dependence 
of the SC order parameter Ai in layer-1 (b), the AF order 
parameter mi in layer-1 (c), and the AF order parameter 
7712 in layer-2 (d) for different gi. The vertical dashed line 
represents Tc in the case of gi = 5.2. The inset in (d) is the 
enlarged drawing of 7712 at low temperatures for gi — 5.0, 5.2. 



Now we examine the case (iii) (Fig. EJa)). The tem- 
perature dependence of the order parameters is shown in 
Figs. [Sl^b)-(d). In this case, we observe quite different 
behaviors of the order parameters compared to the cases 
(i) and (ii). In particular, the superconducting transi- 
tion temperature Tc is always larger than the AF tran- 
sition temperature Tap- There is no coexistence phase 
when Tc < Tap- For Tap < T < Tc, Ai increases as 
T decreases as shown in Fig. [HJb). Below Tap, Ai is 
suppressed. Furthermore, the coexistence phase is lim- 
ited to a finite range of temperature for 6.7 < 51 < 7.1. 



Meanwhile the order parameters mi and m2 increase 
monotonicahy as the temperature decreases as shown in 
Figs. [5jc) and (d). These temperature dependences are 
quahtatively different from the experimentally observed 
one. From this observation one may conclude that it 
is unlikely that there is a coexistence phase of intrinsic 
SC and AF orders in a layer among coupled multilayers. 
Although the coexistence phase of AF and SC orders ap- 
pears in all cases (i)-(iii), the origin of AF order in the 
SC layer is different. What makes the difference between 
the case (iii) and the case (ii) is that in the case (iii) AF 
order in the I = 1 layer with SC is intrinsic order but not 
induced by the other layer. Meanwhile in the case (ii) 
AF order in the I = 1 layer with SC is induced order by 
AF order in the / = 2 layer due to the proximity effect. 

So far we have studied the system with the electron 
hopping restricted to the nearest neighbors. In order to 
examine the system with a realistic Fermi surface ji^r— 
we consider the model with 
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^fc = — 2i(cos kx + cos ky) — 4t' cos k^ cos k 
— 2i"(cos 2kx + cos 2ky) — jjl. 



■y 



(10) 



Here we choose t' /t = -0.20 and t" /t = 0.10. We take 
the chemical potential as ^ ~ —0.839, which corresponds 
to 10% doping in the normal and non-magnetic state. 
The result for the case (i) above with Fi = and V2 = 0.8 
is shown in Fig. |6l The temperature dependence of the 
order parameters is shown in Figs, [nib)-(d) for different 
values of gi. We observe qualitatively similar behaviors 
of the order parameters to the case (i) shown in Fig. |31 
There is a clear enhancement of ?7i2 below the SC tran- 
sition. Intrinsic AF order 7712 appears in the I — 2 layer 
for T < Taf = 0.23 as shown in Fig. E^d). This intrinsic 
order is enhanced below the SC transition. For example, 
there is the SC transition at T = 0.15 for gi = 5.0 as 
shown in Fig. [6l^b). For T < 0.15, m2 is enhanced com- 
pared with the 51 = case. Meanwhile, mi is suppressed 
as shown in Fig. IH^c). The exceptional case is gi = 4.0. 
The temperature dependence of m2 is similar to the other 
cases but the temperature dependence of ttii is qualita- 
tively different. The value of m,i is enhanced below the 
SC transition. This discrepancy is probably associated 
with the change of the Fermi surface shape. 
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FIG. 5. (a) Schematic view of the system for the case (iii) 
(see the text). The temperature dependence of the SC order 
parameter Ai in layer-1 (b), the AF order parameter mi in 
layer-1 (c), and the AF order parameter m2 in layer-2 (d) for 
different gi. 



SUMMARY 



To summarize, we have studied the proximity effect 
and the possibility of coexistence of AF and SC orders 
in a bilayer system. Our mean field theory suggests that 
the experimentally observed enhancement of AF order 
below T^ is associated with the proximity effect. In con- 
trast, if we assume a coexistence phase in a layer among 
coupled multilayers, the temperature dependence of the 
order parameters is qualitatively different from the ex- 
perimentally observed one. 

We believe that this result is not so much affected by 
the shape of the Fermi surface. As we have shown in 



Fig. O the result for a realistic Fermi surface with finite 
doping is qualitatively the same as that for the half-filling 
case. So we expect qualitatively the same proximity effect 
as long as we neglect the possibility of stabilizing other 
orders, such as a charge-density wave or the so-called tt- 
triplet pairing.— i^i^i The absence of the 7r-triplet pairing 
is the unique property of the half-filling case with t' = 
t" = 0. However, there is no experimental evidence for 
the TT-triplet pairing to the best of our knowledge. 
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FIG. 6. (Color online) (a) Schematic view of the system 
for the case (i) with Vi = and V2 = 0.8. The tempera- 
ture dependence of the SC order parameter Ai in layer-1 (b), 
the AF order parameter mi in layer-1 (c), and the AF order 
parameter 7722 in layer-2 (d) for various gi. 
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